Rockfish Historical Harvest Estimates

This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.

Issue	Howard	Bayes
Time series	1999 - present	1977 - present
Bias in SWHS	Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests.	Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases	Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data.	Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations	Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas.	Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation	Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with.	By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Results

• Model: LB_fit_3pH

• Model Run: Gen4int_indcomp_swhsR_FULL_pHB4pars_re0full_altwt_slepsdwn_thru2024_6e+05_2025-11-09 (Naming convention: modelname_thruYEAR_iters_rundate)

• Proportion of parmaters converged (Rhat <= 1.1): 0.9959

• See end of document for summary of unconverged parameters.

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

And black rockfish releases…

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Total Biomass Removal Estimates

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Model fit

Logbook residuals

Figure 20.- Residuals from logbook harvests.

SWHS residuals

Figure 21.- Residuals from SWHS harvests.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

Figure 23.- Mean percent of harvest by charter anglers.

When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(slope|non-pelagic & non-yellowye) For release estimates

## NULL

Weight Fits

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.420354

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.12 while 0.88 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable	area_n	year	user	mean	median	lower	upper	Rhat	sd	skew	area
p_pelagic	13	1	1	0.289	0.271	0.099	0.578	1.019	0.126	0.143	NSEI
p_pelagic	13	1	2	0.116	0.100	0.032	0.287	1.027	0.068	0.234	NSEI
p_pelagic	13	10	1	0.304	0.289	0.113	0.581	1.011	0.121	0.124	NSEI
p_pelagic	13	10	2	0.123	0.108	0.034	0.287	1.015	0.067	0.224	NSEI
p_pelagic	13	11	2	0.128	0.113	0.036	0.300	1.014	0.070	0.214	NSEI
p_pelagic	13	12	1	0.312	0.298	0.113	0.588	1.010	0.126	0.111	NSEI
p_pelagic	13	12	2	0.129	0.114	0.038	0.306	1.012	0.071	0.212	NSEI
p_pelagic	13	13	1	0.316	0.301	0.115	0.580	1.015	0.124	0.121	NSEI
p_pelagic	13	14	2	0.130	0.117	0.039	0.292	1.012	0.066	0.198	NSEI
p_pelagic	13	16	1	0.317	0.303	0.125	0.581	1.016	0.118	0.118	NSEI
p_pelagic	13	17	1	0.325	0.313	0.136	0.581	1.013	0.116	0.103	NSEI
p_pelagic	13	18	1	0.338	0.329	0.140	0.595	1.010	0.117	0.077	NSEI
p_pelagic	13	19	2	0.149	0.137	0.052	0.314	1.012	0.069	0.173	NSEI
p_pelagic	13	2	1	0.292	0.276	0.101	0.579	1.014	0.124	0.129	NSEI
p_pelagic	13	2	2	0.117	0.102	0.033	0.290	1.011	0.066	0.226	NSEI
p_pelagic	13	21	2	0.169	0.156	0.057	0.353	1.010	0.076	0.170	NSEI
p_pelagic	13	3	2	0.118	0.102	0.032	0.298	1.026	0.069	0.233	NSEI
p_pelagic	13	4	2	0.120	0.104	0.033	0.300	1.016	0.070	0.227	NSEI
p_pelagic	13	5	1	0.292	0.272	0.104	0.585	1.012	0.125	0.160	NSEI
p_pelagic	13	5	2	0.119	0.103	0.033	0.284	1.019	0.067	0.240	NSEI
p_pelagic	13	6	1	0.294	0.280	0.103	0.576	1.015	0.126	0.111	NSEI
p_pelagic	13	6	2	0.120	0.104	0.032	0.298	1.019	0.068	0.234	NSEI
p_pelagic	13	7	1	0.295	0.277	0.103	0.577	1.014	0.126	0.143	NSEI
p_pelagic	13	7	2	0.121	0.106	0.033	0.297	1.014	0.068	0.220	NSEI
p_pelagic	13	8	2	0.121	0.106	0.034	0.289	1.022	0.069	0.218	NSEI
p_pelagic	13	9	1	0.305	0.291	0.107	0.576	1.013	0.125	0.112	NSEI
p_pelagic	13	9	2	0.123	0.109	0.036	0.296	1.013	0.067	0.208	NSEI
p_slope	14	13	2	0.073	0.054	0.009	0.255	1.013	0.064	0.298	NSEO
p_pelagic	16	1	1	0.353	0.342	0.115	0.640	1.059	0.138	0.080	SSEO
p_pelagic	16	1	2	0.330	0.314	0.112	0.614	1.060	0.136	0.117	SSEO
p_pelagic	16	10	1	0.358	0.349	0.134	0.637	1.049	0.132	0.068	SSEO
p_pelagic	16	10	2	0.337	0.323	0.120	0.627	1.039	0.133	0.105	SSEO
p_pelagic	16	11	1	0.358	0.347	0.132	0.641	1.049	0.135	0.082	SSEO
p_pelagic	16	11	2	0.342	0.329	0.122	0.630	1.046	0.138	0.095	SSEO
p_pelagic	16	12	1	0.363	0.356	0.135	0.641	1.046	0.134	0.052	SSEO
p_pelagic	16	12	2	0.337	0.326	0.118	0.613	1.033	0.131	0.084	SSEO
p_pelagic	16	13	1	0.360	0.348	0.135	0.658	1.056	0.137	0.087	SSEO
p_pelagic	16	13	2	0.339	0.327	0.121	0.622	1.048	0.132	0.091	SSEO
p_pelagic	16	14	1	0.352	0.343	0.132	0.617	1.033	0.131	0.069	SSEO
p_pelagic	16	14	2	0.345	0.332	0.131	0.619	1.042	0.131	0.099	SSEO
p_pelagic	16	15	1	0.355	0.344	0.143	0.631	1.029	0.128	0.086	SSEO
p_pelagic	16	15	2	0.347	0.333	0.129	0.616	1.031	0.131	0.107	SSEO
p_pelagic	16	16	1	0.353	0.341	0.139	0.629	1.028	0.129	0.093	SSEO
p_pelagic	16	16	2	0.351	0.336	0.139	0.623	1.025	0.129	0.116	SSEO
p_pelagic	16	17	1	0.415	0.406	0.188	0.677	1.033	0.128	0.070	SSEO
p_pelagic	16	17	2	0.367	0.358	0.145	0.637	1.034	0.128	0.070	SSEO
p_pelagic	16	18	1	0.392	0.385	0.174	0.657	1.028	0.129	0.054	SSEO
p_pelagic	16	18	2	0.368	0.359	0.150	0.633	1.033	0.127	0.071	SSEO
p_pelagic	16	19	1	0.382	0.372	0.169	0.633	1.021	0.124	0.080	SSEO
p_pelagic	16	19	2	0.366	0.359	0.150	0.631	1.024	0.125	0.056	SSEO
p_pelagic	16	2	1	0.352	0.343	0.116	0.643	1.057	0.142	0.063	SSEO
p_pelagic	16	2	2	0.330	0.314	0.110	0.613	1.043	0.134	0.119	SSEO
p_pelagic	16	20	1	0.388	0.380	0.176	0.649	1.016	0.123	0.065	SSEO
p_pelagic	16	20	2	0.368	0.358	0.155	0.636	1.019	0.126	0.080	SSEO
p_pelagic	16	21	2	0.370	0.362	0.159	0.623	1.018	0.123	0.065	SSEO
p_pelagic	16	3	1	0.352	0.344	0.112	0.642	1.055	0.142	0.056	SSEO
p_pelagic	16	3	2	0.333	0.321	0.111	0.625	1.048	0.135	0.089	SSEO
p_dsr	16	30	1	0.420	0.419	0.393	0.450	1.028	0.014	0.070	SSEO
p_dsr	16	31	1	0.402	0.401	0.374	0.433	1.013	0.016	0.064	SSEO
p_slope	16	36	1	0.061	0.062	0.041	0.079	1.036	0.009	-0.106	SSEO
p_slope	16	37	1	0.065	0.066	0.046	0.083	1.012	0.009	-0.110	SSEO
p_pelagic	16	4	1	0.353	0.339	0.122	0.656	1.052	0.143	0.098	SSEO
p_pelagic	16	4	2	0.335	0.323	0.111	0.636	1.047	0.137	0.087	SSEO
p_pelagic	16	42	1	0.833	0.833	0.826	0.840	1.055	0.003	0.000	SSEO
p_yellow	16	42	1	0.407	0.407	0.385	0.429	1.022	0.011	0.000	SSEO
p_dsr	16	42	1	0.445	0.438	0.406	0.503	2.015	0.027	0.260	SSEO
p_slope	16	42	1	0.174	0.183	0.115	0.209	2.222	0.028	-0.323	SSEO
p_pelagic	16	5	1	0.356	0.347	0.121	0.634	1.053	0.138	0.065	SSEO
p_pelagic	16	5	2	0.333	0.321	0.115	0.619	1.056	0.135	0.089	SSEO
p_pelagic	16	6	1	0.351	0.338	0.121	0.645	1.071	0.139	0.093	SSEO
p_pelagic	16	6	2	0.335	0.325	0.113	0.613	1.059	0.135	0.074	SSEO
p_pelagic	16	7	1	0.355	0.349	0.123	0.644	1.039	0.138	0.043	SSEO
p_pelagic	16	7	2	0.337	0.325	0.113	0.617	1.046	0.134	0.089	SSEO
p_pelagic	16	8	1	0.356	0.345	0.126	0.651	1.049	0.140	0.079	SSEO
p_pelagic	16	8	2	0.339	0.327	0.112	0.631	1.039	0.135	0.089	SSEO
p_pelagic	16	9	1	0.358	0.343	0.130	0.666	1.050	0.142	0.106	SSEO
p_pelagic	16	9	2	0.338	0.326	0.116	0.625	1.039	0.136	0.088	SSEO
p_pelagic	3	1	1	0.553	0.560	0.286	0.787	1.014	0.133	-0.053	PWSI
p_pelagic	3	1	2	0.187	0.175	0.061	0.380	1.012	0.084	0.142	PWSI
p_pelagic	3	10	1	0.549	0.555	0.275	0.784	1.013	0.133	-0.045	PWSI
p_pelagic	3	10	2	0.189	0.177	0.065	0.391	1.015	0.085	0.141	PWSI
p_pelagic	3	11	1	0.550	0.554	0.281	0.790	1.015	0.132	-0.030	PWSI
p_pelagic	3	12	1	0.551	0.560	0.282	0.785	1.016	0.131	-0.069	PWSI
p_pelagic	3	13	1	0.552	0.559	0.283	0.785	1.010	0.129	-0.054	PWSI
p_pelagic	3	16	2	0.188	0.175	0.067	0.377	1.018	0.082	0.159	PWSI
p_pelagic	3	17	1	0.547	0.550	0.289	0.784	1.012	0.128	-0.023	PWSI
p_pelagic	3	17	2	0.189	0.175	0.065	0.385	1.010	0.083	0.169	PWSI
p_pelagic	3	2	1	0.548	0.554	0.278	0.783	1.013	0.132	-0.045	PWSI
p_pelagic	3	3	1	0.548	0.554	0.271	0.790	1.017	0.135	-0.045	PWSI
p_pelagic	3	4	1	0.552	0.558	0.286	0.787	1.019	0.133	-0.045	PWSI
p_pelagic	3	4	2	0.191	0.178	0.061	0.404	1.013	0.087	0.149	PWSI
p_pelagic	3	5	1	0.553	0.560	0.287	0.781	1.012	0.131	-0.053	PWSI
p_pelagic	3	5	2	0.191	0.178	0.062	0.392	1.022	0.086	0.151	PWSI
p_pelagic	3	6	1	0.551	0.556	0.281	0.785	1.019	0.129	-0.039	PWSI
p_pelagic	3	6	2	0.189	0.175	0.065	0.395	1.014	0.085	0.165	PWSI
p_pelagic	3	7	1	0.549	0.554	0.282	0.781	1.011	0.132	-0.038	PWSI
p_pelagic	3	7	2	0.188	0.176	0.059	0.380	1.013	0.083	0.144	PWSI
p_pelagic	3	8	1	0.551	0.555	0.295	0.788	1.013	0.132	-0.030	PWSI
p_pelagic	3	8	2	0.188	0.176	0.062	0.391	1.010	0.085	0.141	PWSI
p_pelagic	3	9	1	0.553	0.558	0.276	0.796	1.018	0.134	-0.037	PWSI
p_pelagic	3	9	2	0.189	0.175	0.064	0.383	1.015	0.083	0.168	PWSI
p_pelagic	4	1	1	0.546	0.551	0.284	0.793	1.012	0.132	-0.038	PWSO
p_pelagic	4	10	1	0.549	0.552	0.292	0.789	1.013	0.130	-0.023	PWSO
p_pelagic	4	10	2	0.554	0.558	0.294	0.803	1.015	0.132	-0.030	PWSO
p_pelagic	4	11	1	0.550	0.557	0.290	0.780	1.015	0.130	-0.054	PWSO
p_pelagic	4	11	2	0.558	0.560	0.285	0.795	1.010	0.132	-0.015	PWSO
p_pelagic	4	12	1	0.546	0.549	0.288	0.781	1.015	0.129	-0.023	PWSO
p_pelagic	4	12	2	0.554	0.558	0.287	0.798	1.014	0.131	-0.031	PWSO
p_pelagic	4	13	1	0.545	0.551	0.283	0.789	1.013	0.131	-0.046	PWSO
p_pelagic	4	13	2	0.552	0.558	0.279	0.797	1.019	0.132	-0.045	PWSO
p_pelagic	4	15	2	0.551	0.556	0.293	0.784	1.015	0.129	-0.039	PWSO
p_pelagic	4	16	1	0.544	0.551	0.284	0.771	1.020	0.130	-0.054	PWSO
p_pelagic	4	16	2	0.556	0.559	0.303	0.793	1.017	0.132	-0.023	PWSO
p_pelagic	4	17	2	0.550	0.553	0.287	0.794	1.020	0.131	-0.023	PWSO
p_pelagic	4	18	1	0.546	0.552	0.296	0.789	1.014	0.130	-0.046	PWSO
p_pelagic	4	18	2	0.556	0.563	0.292	0.800	1.017	0.134	-0.052	PWSO
p_pelagic	4	19	2	0.551	0.551	0.295	0.791	1.012	0.131	0.000	PWSO
p_pelagic	4	2	1	0.549	0.552	0.286	0.785	1.014	0.130	-0.023	PWSO
p_pelagic	4	3	1	0.552	0.557	0.281	0.793	1.014	0.132	-0.038	PWSO
p_pelagic	4	3	2	0.552	0.556	0.288	0.796	1.020	0.133	-0.030	PWSO
p_pelagic	4	4	2	0.551	0.556	0.288	0.793	1.019	0.131	-0.038	PWSO
p_pelagic	4	5	1	0.543	0.547	0.285	0.783	1.016	0.132	-0.030	PWSO
p_pelagic	4	5	2	0.555	0.560	0.296	0.795	1.019	0.133	-0.038	PWSO
p_pelagic	4	6	1	0.552	0.560	0.290	0.790	1.010	0.133	-0.060	PWSO
p_pelagic	4	7	1	0.549	0.553	0.288	0.791	1.012	0.131	-0.031	PWSO
p_pelagic	4	7	2	0.555	0.559	0.293	0.790	1.018	0.131	-0.030	PWSO
p_pelagic	4	8	1	0.551	0.555	0.280	0.788	1.011	0.131	-0.031	PWSO
p_pelagic	4	8	2	0.552	0.557	0.280	0.791	1.013	0.135	-0.037	PWSO
p_pelagic	4	9	1	0.551	0.556	0.288	0.795	1.018	0.134	-0.037	PWSO
p_pelagic	4	9	2	0.552	0.556	0.290	0.797	1.013	0.133	-0.030	PWSO
p_black	8	40	1	0.475	0.415	0.341	0.963	1.016	0.178	0.337	afognak

Table Xb. True parameters with an Rhat in excess of 1.1

variable	area_n	year/species	mean	median	lower	upper	Rhat	sd	skew	year	area
beta2_pH	10	1	1.012	0.670	0.287	4.242	1.193	1.047	0.327		northeast
beta0_pH	13	4	0.603	0.692	-0.964	1.676	1.141	0.664	-0.134		NSEI
beta1_pH	13	4	1.396	1.341	0.234	2.904	1.138	0.677	0.081		NSEI
beta0_pH	14	4	0.399	0.453	-0.762	1.286	1.129	0.530	-0.102		NSEO
beta1_pH	14	4	1.967	1.921	0.911	3.246	1.126	0.603	0.076		NSEO
beta0_pH	15	4	0.573	0.692	-1.026	1.548	1.314	0.650	-0.183		SSEI
beta1_pH	15	4	1.405	1.296	0.314	3.006	1.291	0.691	0.158		SSEI
beta0_yellow	15		-1.282	-1.237	-1.967	-1.027	1.254	0.259	-0.174		SSEI
beta1_yellow	15		0.192	0.073	0.011	1.081	1.156	0.297	0.401		SSEI
beta1_pH	16	4	3.529	3.520	2.216	4.964	1.107	0.699	0.013		SSEO
beta0_pelagic	16		-0.689	-0.627	-1.624	-0.013	1.161	0.454	-0.136		SSEO
beta1_pelagic	16		1.890	1.865	0.832	3.055	1.181	0.634	0.039		SSEO
mu_beta2_pH	2	1	0.946	0.588	0.061	3.962	1.115	1.006	0.356		NG
beta2_yellow	2		-5.744	-3.671	-20.442	-0.264	1.338	5.590	-0.371		NG
mu_beta0_pH	3	4	0.690	0.757	-0.543	1.622	1.136	0.559	-0.120		PWSI
tau_beta1_pH	3	4	5.778	2.866	0.231	34.183	1.137	9.370	0.311		PWSI
beta1_pelagic	3		0.769	0.718	0.309	1.587	1.108	0.301	0.169		PWSI
beta0_pelagic	4		0.219	0.231	-0.353	0.668	1.107	0.243	-0.049		PWSO
beta1_pelagic	4		0.832	0.822	0.334	1.506	1.112	0.268	0.037		PWSO
beta2_pH	8	1	0.683	0.344	0.123	4.197	1.132	1.045	0.324		afognak
beta2_pH	9	1	0.899	0.512	0.148	4.022	1.127	1.037	0.373		eastside
beta2_pH	9	3	0.282	0.195	0.088	0.865	1.125	0.520	0.167		eastside
lambda_H	15		0.019	0.014	0.003	0.061	1.226	0.043	0.117		SSEI

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic	non.BSAI.SOKO.Harv	BSAI.SOKO.Harv	non.BSAI.SOKO.Rel	BSAI.SOKO.Rel
n unconverged at Rhat = 1.01	449.00	200.00	1,016.00	342.00
Prop. of total	0.69	0.31	0.75	0.25
Prop. median less than 1	0.27	0.26	0.35	0.45
Prop. median less than 10	0.42	0.56	0.58	0.76
Prop. median less than 50	0.60	0.81	0.75	0.85
Median of the median values	24.67	6.32	4.22	1.48
Mean of the median values	348.75	134.02	184.73	105.38
Mean skew of posteriors	0.26	0.28	0.27	0.24

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic	non.BSAI.SOKO.Harv	BSAI.SOKO.Harv	non.BSAI.SOKO.Rel	BSAI.SOKO.Rel
n unconverged at Rhat = 1.1	12.00	14.00	62.00	47.00
Prop. of total	0.46	0.54	0.57	0.43
Prop. median less than 1	0.42	0.43	0.71	0.66
Prop. median less than 10	0.42	1.00	0.84	0.89
Prop. median less than 50	0.42	1.00	0.87	0.96
Median of the median values	180.53	1.02	0.06	0.33
Mean of the median values	421.52	2.29	37.09	4.61
Mean skew of posteriors	0.23	0.19	0.16	0.17

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable	Freq
beta0_pelagic	1
beta0_pH	3
beta0_yellow	1
beta1_pelagic	2
beta1_pH	3
beta1_yellow	1
beta2_pH	4
beta2_yellow	1
eps_dsr	1
eps_pH	14
eps_sl	1
H_ay	4
H_ayg	1
H_ayu	4
Hb_ay	4
Hb_ayg	1
Hb_ayu	4
Hd_ayg	1
Ho_ayg	1
Ho_ayu	3
Hp_ay	4
Hp_ayg	1
Hp_ayu	4
Hs_ay	1
Hs_ayg	1
Hs_ayu	2
Hy_ay	4
Hy_ayg	1
Hy_ayu	8
lambda_H	1
mean_eps_pH	1
mu_beta0_pH	1
mu_beta2_pH	1
p_dsr	1
p_slope	1
pH	9
R_ay	4
R_ayg	2
R_ayu	5
Rb_ay	4
Rb_ayg	2
Rb_ayu	5
Rd_ay	2
Rd_ayg	3
Rd_ayu	3
re_dsr	1
re_pH	10
re_slope	1
Ro_ay	24
Ro_ayg	13
Ro_ayu	39
Rp_ay	4
Rp_ayg	2
Rp_ayu	4
Rs_ay	3
Rs_ayg	10
Rs_ayu	20
Ry_ay	10
Ry_ayg	3
Ry_ayu	25
tau_beta1_pH	1

Parameter estimates:

Summary Table of Parameter Estimates

Parameter	mean	sd	Lower_CI	Median	Upper_CI
mu_bc_H[1]	-0.020	0.072	-0.149	-0.024	0.129
mu_bc_H[2]	0.075	0.053	-0.022	0.072	0.186
mu_bc_H[3]	-0.266	0.068	-0.394	-0.268	-0.129
mu_bc_H[4]	-0.931	0.183	-1.312	-0.926	-0.580
mu_bc_H[5]	1.355	1.197	0.129	1.082	4.330
mu_bc_H[6]	-2.375	0.321	-2.994	-2.381	-1.750
mu_bc_H[7]	-0.453	0.103	-0.664	-0.449	-0.261
mu_bc_H[8]	0.245	0.350	-0.334	0.207	1.023
mu_bc_H[9]	-0.247	0.124	-0.487	-0.251	0.002
mu_bc_H[10]	0.060	0.065	-0.065	0.059	0.190
mu_bc_H[11]	-0.128	0.037	-0.202	-0.128	-0.060
mu_bc_H[12]	-0.256	0.111	-0.491	-0.251	-0.045
mu_bc_H[13]	-0.162	0.072	-0.297	-0.163	-0.023
mu_bc_H[14]	-0.306	0.092	-0.483	-0.304	-0.124
mu_bc_H[15]	-0.336	0.046	-0.427	-0.337	-0.243
mu_bc_H[16]	-0.176	0.366	-0.828	-0.202	0.574
mu_bc_R[1]	1.036	0.123	0.793	1.034	1.272
mu_bc_R[2]	1.011	0.075	0.859	1.013	1.149
mu_bc_R[3]	0.227	0.093	0.038	0.229	0.406
mu_bc_R[4]	0.662	0.199	0.261	0.670	1.024
mu_bc_R[5]	1.674	0.374	0.965	1.673	2.452
mu_bc_R[6]	-1.412	0.338	-2.067	-1.416	-0.737
mu_bc_R[7]	0.511	0.170	0.171	0.515	0.837
mu_bc_R[8]	0.647	0.144	0.359	0.645	0.936
mu_bc_R[9]	0.385	0.171	0.023	0.394	0.698
mu_bc_R[10]	0.897	0.100	0.707	0.896	1.092
mu_bc_R[11]	1.110	0.074	0.967	1.109	1.262
mu_bc_R[12]	0.867	0.164	0.539	0.869	1.188
mu_bc_R[13]	1.032	0.091	0.849	1.033	1.208
mu_bc_R[14]	0.862	0.123	0.620	0.863	1.116
mu_bc_R[15]	0.804	0.097	0.614	0.804	0.992
mu_bc_R[16]	1.353	0.114	1.118	1.354	1.571
tau_pH[1]	2.453	0.175	2.124	2.452	2.806
tau_pH[2]	1.607	0.151	1.332	1.600	1.919
tau_pH[3]	1.193	0.108	0.992	1.189	1.411
tau_pH[4]	1.026	0.147	0.771	1.011	1.342
tau_pH[5]	0.422	0.075	0.295	0.414	0.583
beta0_pH[1,1]	0.551	0.166	0.212	0.556	0.868
beta0_pH[2,1]	1.428	0.164	1.063	1.438	1.727
beta0_pH[3,1]	1.410	0.208	0.928	1.434	1.765
beta0_pH[4,1]	1.595	0.219	1.125	1.605	1.998
beta0_pH[5,1]	-0.819	0.407	-1.824	-0.777	-0.134
beta0_pH[6,1]	-0.588	0.521	-2.107	-0.445	0.043
beta0_pH[7,1]	0.757	0.010	0.736	0.757	0.777
beta0_pH[8,1]	-0.652	0.410	-1.776	-0.581	-0.092
beta0_pH[9,1]	-0.341	0.312	-1.020	-0.327	0.208
beta0_pH[10,1]	0.245	0.205	-0.190	0.253	0.599
beta0_pH[11,1]	0.147	0.304	-0.548	0.171	0.662
beta0_pH[12,1]	0.627	0.186	0.250	0.632	0.973
beta0_pH[13,1]	-0.657	0.144	-0.964	-0.648	-0.397
beta0_pH[14,1]	0.014	0.175	-0.339	0.015	0.352
beta0_pH[15,1]	0.072	0.175	-0.307	0.086	0.381
beta0_pH[16,1]	-0.854	0.214	-1.265	-0.863	-0.435
beta0_pH[1,2]	2.564	0.085	2.407	2.564	2.737
beta0_pH[2,2]	3.079	0.022	3.038	3.079	3.121
beta0_pH[3,2]	3.626	0.041	3.548	3.625	3.706
beta0_pH[4,2]	2.936	0.019	2.899	2.936	2.975
beta0_pH[5,2]	4.569	1.529	2.745	4.252	8.360
beta0_pH[6,2]	3.324	0.157	3.020	3.318	3.633
beta0_pH[7,2]	2.032	0.051	1.934	2.031	2.131
beta0_pH[8,2]	2.824	0.065	2.700	2.823	2.954
beta0_pH[9,2]	3.954	0.210	3.580	3.945	4.397
beta0_pH[10,2]	3.997	0.127	3.753	3.994	4.260
beta0_pH[11,2]	1.885	0.009	1.867	1.885	1.903
beta0_pH[12,2]	1.792	0.068	1.664	1.788	1.927
beta0_pH[13,2]	2.496	0.026	2.445	2.496	2.547
beta0_pH[14,2]	1.803	0.024	1.756	1.803	1.850
beta0_pH[15,2]	2.545	0.017	2.513	2.545	2.578
beta0_pH[16,2]	2.766	0.025	2.718	2.766	2.816
beta0_pH[1,3]	1.071	0.395	0.180	1.143	1.661
beta0_pH[2,3]	2.126	0.017	2.091	2.127	2.160
beta0_pH[3,3]	2.692	0.024	2.644	2.692	2.739
beta0_pH[4,3]	2.981	0.035	2.913	2.980	3.046
beta0_pH[5,3]	3.502	2.411	0.897	2.921	10.112
beta0_pH[6,3]	1.363	0.376	0.731	1.328	2.183
beta0_pH[7,3]	0.859	0.050	0.760	0.858	0.956
beta0_pH[8,3]	0.179	0.611	-1.509	0.290	0.631
beta0_pH[9,3]	-2.299	1.038	-4.636	-2.153	-0.577
beta0_pH[10,3]	-0.497	0.834	-2.386	-0.411	0.861
beta0_pH[11,3]	1.088	0.550	-0.139	1.125	2.037
beta0_pH[12,3]	0.963	0.652	-0.466	0.975	2.169
beta0_pH[13,3]	0.983	0.533	-0.249	1.021	1.941
beta0_pH[14,3]	0.563	0.464	-0.482	0.609	1.357
beta0_pH[15,3]	0.722	0.443	-0.240	0.741	1.500
beta0_pH[16,3]	1.178	0.634	-0.053	1.170	2.440
beta0_pH[11,4]	0.859	0.548	-0.408	0.908	1.826
beta0_pH[12,4]	0.905	0.692	-0.503	0.923	2.220
beta0_pH[13,4]	0.603	0.664	-0.964	0.692	1.676
beta0_pH[14,4]	0.399	0.530	-0.762	0.453	1.286
beta0_pH[15,4]	0.573	0.650	-1.026	0.692	1.548
beta0_pH[16,4]	0.862	0.621	-0.477	0.892	2.026
beta0_pH[11,5]	2.964	0.145	2.685	2.964	3.258
beta0_pH[12,5]	-0.919	0.509	-1.966	-0.902	0.025
beta0_pH[13,5]	4.922	0.426	4.165	4.906	5.832
beta0_pH[14,5]	3.911	0.374	3.226	3.888	4.703
beta0_pH[15,5]	2.289	0.129	2.052	2.286	2.554
beta0_pH[16,5]	0.244	0.063	0.158	0.232	0.401
beta1_pH[1,1]	2.573	0.448	1.819	2.533	3.589
beta1_pH[2,1]	1.577	0.348	1.035	1.537	2.535
beta1_pH[3,1]	1.782	0.391	1.135	1.743	2.704
beta1_pH[4,1]	1.870	0.407	1.197	1.826	2.863
beta1_pH[5,1]	2.359	0.523	1.517	2.299	3.671
beta1_pH[6,1]	2.270	0.875	1.111	2.039	4.440
beta1_pH[7,1]	0.000	0.000	0.000	0.000	0.000
beta1_pH[8,1]	2.451	0.871	1.217	2.240	4.469
beta1_pH[9,1]	2.121	0.698	1.160	1.985	4.245
beta1_pH[10,1]	2.213	0.372	1.575	2.190	3.035
beta1_pH[11,1]	3.266	0.518	2.447	3.179	4.489
beta1_pH[12,1]	2.900	0.274	2.369	2.902	3.454
beta1_pH[13,1]	3.457	0.237	3.014	3.448	3.955
beta1_pH[14,1]	2.834	0.321	2.227	2.827	3.496
beta1_pH[15,1]	2.383	0.260	1.901	2.369	2.920
beta1_pH[16,1]	4.004	0.297	3.417	4.023	4.562
beta1_pH[1,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[2,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[3,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[4,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[5,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[6,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[7,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[8,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[9,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[10,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[11,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[12,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[13,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[14,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[15,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[16,2]	0.000	0.000	0.000	0.000	0.000
beta1_pH[1,3]	1.950	1.246	0.064	1.780	5.349
beta1_pH[2,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[3,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[4,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[5,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[6,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[7,3]	0.000	0.000	0.000	0.000	0.000
beta1_pH[8,3]	4.773	2.577	0.790	4.279	11.532
beta1_pH[9,3]	7.651	2.637	3.637	7.243	13.594
beta1_pH[10,3]	6.589	2.670	2.732	6.127	12.618
beta1_pH[11,3]	2.201	0.597	1.120	2.165	3.493
beta1_pH[12,3]	3.250	0.674	2.031	3.225	4.741
beta1_pH[13,3]	1.992	0.566	0.920	1.971	3.293
beta1_pH[14,3]	2.210	0.518	1.298	2.178	3.364
beta1_pH[15,3]	2.102	0.513	1.164	2.082	3.190
beta1_pH[16,3]	3.705	0.689	2.291	3.715	5.009
beta1_pH[11,4]	1.923	0.593	0.820	1.904	3.193
beta1_pH[12,4]	3.178	0.711	1.805	3.157	4.618
beta1_pH[13,4]	1.396	0.677	0.234	1.341	2.904
beta1_pH[14,4]	1.967	0.603	0.911	1.921	3.246
beta1_pH[15,4]	1.405	0.691	0.314	1.296	3.006
beta1_pH[16,4]	3.529	0.699	2.216	3.520	4.964
beta1_pH[11,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[12,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[13,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[14,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[15,5]	0.000	0.000	0.000	0.000	0.000
beta1_pH[16,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[1,1]	1.255	0.851	0.360	0.931	3.384
beta2_pH[2,1]	1.631	0.968	0.239	1.473	3.427
beta2_pH[3,1]	1.412	0.837	0.229	1.270	2.941
beta2_pH[4,1]	1.184	0.790	0.229	0.941	2.917
beta2_pH[5,1]	1.879	1.612	0.260	1.425	6.419
beta2_pH[6,1]	0.879	1.375	0.106	0.367	5.103
beta2_pH[7,1]	0.000	0.000	0.000	0.000	0.000
beta2_pH[8,1]	0.683	1.045	0.123	0.344	4.197
beta2_pH[9,1]	0.899	1.037	0.148	0.512	4.022
beta2_pH[10,1]	1.012	1.047	0.287	0.670	4.242
beta2_pH[11,1]	0.463	0.264	0.188	0.400	1.104
beta2_pH[12,1]	1.556	0.792	0.647	1.301	3.632
beta2_pH[13,1]	1.002	0.388	0.550	0.923	2.065
beta2_pH[14,1]	1.144	0.671	0.435	0.945	3.234
beta2_pH[15,1]	1.052	0.632	0.374	0.877	3.024
beta2_pH[16,1]	1.052	0.310	0.644	0.995	1.773
beta2_pH[1,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[2,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[3,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[4,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[5,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[6,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[7,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[8,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[9,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[10,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[11,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[12,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[13,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[14,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[15,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[16,2]	0.000	0.000	0.000	0.000	0.000
beta2_pH[1,3]	6.870	6.654	0.173	5.144	22.973
beta2_pH[2,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[3,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[4,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[5,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[6,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[7,3]	0.000	0.000	0.000	0.000	0.000
beta2_pH[8,3]	1.449	1.367	0.070	1.020	5.630
beta2_pH[9,3]	0.282	0.520	0.088	0.195	0.865
beta2_pH[10,3]	0.406	0.604	0.105	0.260	1.699
beta2_pH[11,3]	-6.178	4.915	-19.394	-4.832	-0.583
beta2_pH[12,3]	-6.983	4.798	-19.422	-5.802	-1.229
beta2_pH[13,3]	-6.297	4.833	-19.064	-4.843	-0.959
beta2_pH[14,3]	-6.487	4.860	-20.336	-5.024	-1.217
beta2_pH[15,3]	-6.575	4.828	-19.335	-5.232	-1.100
beta2_pH[16,3]	-6.495	4.745	-19.410	-5.048	-1.479
beta2_pH[11,4]	-5.876	4.424	-17.149	-4.717	-0.776
beta2_pH[12,4]	-6.192	4.367	-16.798	-5.033	-1.059
beta2_pH[13,4]	-6.010	4.393	-16.858	-4.762	-0.901
beta2_pH[14,4]	-6.081	4.249	-16.784	-4.973	-1.074
beta2_pH[15,4]	-5.845	4.326	-15.949	-4.741	-0.802
beta2_pH[16,4]	-6.169	4.126	-16.574	-4.996	-1.396
beta2_pH[11,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[12,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[13,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[14,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[15,5]	0.000	0.000	0.000	0.000	0.000
beta2_pH[16,5]	0.000	0.000	0.000	0.000	0.000
beta3_pH[1,1]	34.940	1.110	32.887	34.923	37.356
beta3_pH[2,1]	32.630	1.471	30.094	32.514	36.108
beta3_pH[3,1]	33.069	1.286	30.524	33.022	35.758
beta3_pH[4,1]	32.730	1.763	29.737	32.639	36.619
beta3_pH[5,1]	27.643	1.260	25.948	27.473	31.207
beta3_pH[6,1]	33.229	3.321	26.351	33.062	40.720
beta3_pH[7,1]	25.000	0.000	25.000	25.000	25.000
beta3_pH[8,1]	35.151	3.803	27.261	35.274	42.105
beta3_pH[9,1]	32.403	2.827	28.073	31.939	39.765
beta3_pH[10,1]	32.636	1.271	30.330	32.587	35.336
beta3_pH[11,1]	31.566	1.110	29.453	31.566	33.944
beta3_pH[12,1]	30.787	0.483	29.784	30.810	31.728
beta3_pH[13,1]	31.415	0.548	30.361	31.398	32.489
beta3_pH[14,1]	32.483	0.782	31.188	32.411	34.148
beta3_pH[15,1]	31.013	0.727	29.521	31.015	32.439
beta3_pH[16,1]	30.145	0.345	29.443	30.141	30.849
beta3_pH[1,2]	34.608	5.498	25.405	34.809	43.471
beta3_pH[2,2]	34.507	5.560	25.468	34.469	43.540
beta3_pH[3,2]	30.573	7.150	18.644	30.771	42.310
beta3_pH[4,2]	30.524	7.141	18.755	30.697	42.315
beta3_pH[5,2]	30.499	7.185	18.622	30.687	42.430
beta3_pH[6,2]	30.462	7.226	18.648	30.435	42.348
beta3_pH[7,2]	25.000	0.000	25.000	25.000	25.000
beta3_pH[8,2]	25.000	0.000	25.000	25.000	25.000
beta3_pH[9,2]	25.000	0.000	25.000	25.000	25.000
beta3_pH[10,2]	25.000	0.000	25.000	25.000	25.000
beta3_pH[11,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[12,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[13,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[14,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[15,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[16,2]	24.000	0.000	24.000	24.000	24.000
beta3_pH[1,3]	38.072	3.398	29.390	39.357	42.795
beta3_pH[2,3]	34.487	5.421	25.412	34.473	43.573
beta3_pH[3,3]	30.429	7.129	18.738	30.273	42.360
beta3_pH[4,3]	30.641	7.205	18.600	30.820	42.284
beta3_pH[5,3]	25.000	0.000	25.000	25.000	25.000
beta3_pH[6,3]	25.000	0.000	25.000	25.000	25.000
beta3_pH[7,3]	25.000	0.000	25.000	25.000	25.000
beta3_pH[8,3]	41.959	1.467	37.575	42.056	43.886
beta3_pH[9,3]	36.184	1.900	34.525	35.437	41.476
beta3_pH[10,3]	37.979	2.548	34.668	37.339	43.474
beta3_pH[11,3]	40.750	1.050	38.895	40.772	42.705
beta3_pH[12,3]	42.039	0.712	40.322	42.285	42.832
beta3_pH[13,3]	41.470	0.752	39.969	41.481	42.876
beta3_pH[14,3]	40.515	0.697	39.083	40.554	41.786
beta3_pH[15,3]	40.646	0.642	39.199	40.670	41.886
beta3_pH[16,3]	41.840	0.392	41.034	41.897	42.608
beta3_pH[11,4]	40.655	1.112	38.478	40.581	42.994
beta3_pH[12,4]	41.806	0.819	40.228	42.080	42.932
beta3_pH[13,4]	41.881	1.233	39.037	41.984	43.807
beta3_pH[14,4]	40.440	0.767	38.840	40.470	41.839
beta3_pH[15,4]	41.405	1.204	38.620	41.398	43.605
beta3_pH[16,4]	41.655	0.428	40.906	41.663	42.559
beta3_pH[11,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[12,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[13,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[14,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[15,5]	25.000	0.000	25.000	25.000	25.000
beta3_pH[16,5]	25.000	0.000	25.000	25.000	25.000
beta4_pH[1,1]	0.264	0.020	0.224	0.264	0.302
beta4_pH[2,1]	0.211	0.005	0.201	0.211	0.221
beta4_pH[3,1]	0.173	0.007	0.159	0.173	0.187
beta4_pH[4,1]	0.293	0.011	0.272	0.294	0.315
beta4_pH[5,1]	0.426	0.266	0.091	0.386	0.948
beta4_pH[6,1]	0.779	0.368	0.310	0.729	1.480
beta4_pH[7,1]	0.855	0.348	0.331	0.808	1.671
beta4_pH[8,1]	0.157	0.069	0.044	0.150	0.315
beta4_pH[9,1]	0.213	0.086	0.060	0.209	0.390
beta4_pH[10,1]	0.001	0.001	0.000	0.001	0.003
beta4_pH[11,1]	0.447	0.012	0.424	0.447	0.469
beta4_pH[12,1]	0.642	0.022	0.599	0.641	0.686
beta4_pH[13,1]	0.631	0.013	0.605	0.630	0.657
beta4_pH[14,1]	0.670	0.032	0.610	0.669	0.735
beta4_pH[15,1]	0.798	0.021	0.758	0.797	0.838
beta4_pH[16,1]	1.220	0.043	1.141	1.220	1.307
beta4_pH[1,2]	0.309	0.096	0.137	0.304	0.511
beta4_pH[2,2]	0.515	0.023	0.472	0.514	0.561
beta4_pH[3,2]	0.446	0.012	0.423	0.445	0.468
beta4_pH[4,2]	0.654	0.038	0.582	0.654	0.730
beta4_pH[5,2]	0.738	0.953	0.112	0.535	2.746
beta4_pH[6,2]	0.928	0.828	0.210	0.745	2.762
beta4_pH[7,2]	0.900	0.445	0.335	0.818	1.968
beta4_pH[8,2]	0.843	0.201	0.480	0.833	1.269
beta4_pH[9,2]	1.161	1.053	0.302	0.891	3.940
beta4_pH[10,2]	0.455	0.052	0.356	0.454	0.560
beta4_pH[11,2]	1.267	0.050	1.171	1.268	1.364
beta4_pH[12,2]	0.584	0.110	0.390	0.579	0.806
beta4_pH[13,2]	0.945	0.044	0.857	0.945	1.033
beta4_pH[14,2]	0.931	0.062	0.817	0.930	1.054
beta4_pH[15,2]	1.011	0.026	0.961	1.011	1.062
beta4_pH[16,2]	1.406	0.075	1.265	1.405	1.551
beta4_pH[1,3]	0.284	0.086	0.131	0.281	0.464
beta4_pH[2,3]	0.051	0.017	0.019	0.050	0.088
beta4_pH[3,3]	0.052	0.010	0.034	0.052	0.071
beta4_pH[4,3]	0.173	0.032	0.113	0.173	0.235
beta4_pH[5,3]	0.747	0.774	0.123	0.570	2.473
beta4_pH[6,3]	0.885	0.827	0.212	0.717	2.549
beta4_pH[7,3]	0.900	0.463	0.343	0.811	1.956
beta4_pH[8,3]	0.415	0.200	0.110	0.388	0.867
beta4_pH[9,3]	0.930	0.908	0.158	0.708	3.078
beta4_pH[10,3]	0.061	0.042	0.004	0.056	0.158
beta4_pH[11,3]	0.364	0.016	0.333	0.364	0.397
beta4_pH[12,3]	0.440	0.041	0.363	0.438	0.525
beta4_pH[13,3]	0.840	0.037	0.767	0.840	0.912
beta4_pH[14,3]	0.825	0.068	0.700	0.825	0.962
beta4_pH[15,3]	0.635	0.017	0.603	0.635	0.668
beta4_pH[16,3]	0.638	0.036	0.569	0.638	0.710
beta4_pH[1,4]	0.630	0.336	0.220	0.559	1.442
beta4_pH[2,4]	0.803	0.958	0.119	0.580	2.902
beta4_pH[3,4]	0.835	0.861	0.109	0.587	2.993
beta4_pH[4,4]	0.730	0.487	0.218	0.624	1.909
beta4_pH[5,4]	0.775	0.924	0.122	0.578	2.791
beta4_pH[6,4]	0.902	0.785	0.210	0.728	2.618
beta4_pH[7,4]	0.897	0.409	0.350	0.827	1.879
beta4_pH[8,4]	0.753	0.611	0.154	0.600	2.232
beta4_pH[9,4]	0.770	0.807	0.129	0.567	2.650
beta4_pH[10,4]	1.470	3.035	0.008	0.394	10.910
beta4_pH[11,4]	0.444	0.017	0.411	0.444	0.478
beta4_pH[12,4]	0.561	0.055	0.457	0.560	0.674
beta4_pH[13,4]	0.925	0.036	0.855	0.926	0.995
beta4_pH[14,4]	0.922	0.074	0.778	0.922	1.071
beta4_pH[15,4]	0.714	0.020	0.675	0.714	0.755
beta4_pH[16,4]	0.689	0.042	0.610	0.689	0.776
beta4_pH[1,5]	0.628	0.343	0.209	0.566	1.419
beta4_pH[2,5]	0.795	0.801	0.116	0.574	2.852
beta4_pH[3,5]	0.856	0.931	0.112	0.596	3.187
beta4_pH[4,5]	0.727	0.529	0.213	0.622	1.797
beta4_pH[5,5]	0.770	0.826	0.114	0.561	2.755
beta4_pH[6,5]	0.907	0.826	0.210	0.723	2.758
beta4_pH[7,5]	0.889	0.433	0.348	0.814	1.895
beta4_pH[8,5]	0.763	0.864	0.157	0.589	2.402
beta4_pH[9,5]	0.790	0.984	0.120	0.573	2.699
beta4_pH[10,5]	1.488	3.110	0.007	0.402	10.875
beta4_pH[11,5]	0.602	0.052	0.502	0.601	0.707
beta4_pH[12,5]	0.580	0.164	0.308	0.567	0.955
beta4_pH[13,5]	0.773	0.063	0.655	0.771	0.907
beta4_pH[14,5]	0.846	0.119	0.635	0.838	1.106
beta4_pH[15,5]	1.144	0.076	1.005	1.141	1.295
beta4_pH[16,5]	19.106	3.875	11.029	19.646	24.729
beta0_pelagic[1]	2.267	0.036	2.200	2.266	2.338
beta0_pelagic[2]	1.538	0.007	1.524	1.538	1.552
beta0_pelagic[3]	0.186	0.259	-0.446	0.244	0.563
beta0_pelagic[4]	0.219	0.243	-0.353	0.231	0.668
beta0_pelagic[5]	1.448	0.154	1.158	1.445	1.764
beta0_pelagic[6]	1.353	0.060	1.240	1.351	1.475
beta0_pelagic[7]	1.480	0.039	1.407	1.479	1.557
beta0_pelagic[8]	1.673	0.044	1.593	1.673	1.762
beta0_pelagic[9]	2.794	0.066	2.667	2.793	2.928
beta0_pelagic[10]	2.603	0.031	2.548	2.601	2.666
beta0_pelagic[11]	-0.098	0.141	-0.400	-0.092	0.164
beta0_pelagic[12]	1.801	0.035	1.738	1.799	1.873
beta0_pelagic[13]	-1.061	0.416	-1.711	-1.096	-0.202
beta0_pelagic[14]	-0.280	0.168	-0.644	-0.269	0.018
beta0_pelagic[15]	-0.211	0.034	-0.298	-0.197	-0.180
beta0_pelagic[16]	-0.689	0.454	-1.624	-0.627	-0.013
beta1_pelagic[1]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[2]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[3]	0.769	0.301	0.309	0.718	1.587
beta1_pelagic[4]	0.832	0.268	0.334	0.822	1.506
beta1_pelagic[5]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[6]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[7]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[8]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[9]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[10]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[11]	1.104	0.220	0.731	1.082	1.641
beta1_pelagic[12]	0.000	0.000	0.000	0.000	0.000
beta1_pelagic[13]	2.885	0.608	1.648	2.860	3.915
beta1_pelagic[14]	1.758	0.434	1.091	1.671	2.823
beta1_pelagic[15]	0.094	0.161	0.000	0.016	0.609
beta1_pelagic[16]	1.890	0.634	0.832	1.865	3.055
beta2_pelagic[1]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[2]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[3]	6.139	5.921	0.169	4.433	20.583
beta2_pelagic[4]	7.403	6.177	0.317	5.985	22.792
beta2_pelagic[5]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[6]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[7]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[8]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[9]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[10]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[11]	1.881	1.526	0.284	1.389	5.920
beta2_pelagic[12]	0.000	0.000	0.000	0.000	0.000
beta2_pelagic[13]	0.175	0.070	0.095	0.159	0.358
beta2_pelagic[14]	0.924	0.950	0.217	0.616	3.624
beta2_pelagic[15]	2.665	1.901	0.165	2.308	7.240
beta2_pelagic[16]	0.534	0.793	0.112	0.261	3.174
beta3_pelagic[1]	28.851	3.860	22.298	28.939	35.120
beta3_pelagic[2]	28.888	3.869	22.290	28.935	35.145
beta3_pelagic[3]	27.807	2.664	23.202	28.094	32.028
beta3_pelagic[4]	26.163	1.529	23.060	26.280	29.848
beta3_pelagic[5]	29.449	6.552	18.608	29.325	40.396
beta3_pelagic[6]	29.476	6.606	18.627	29.473	40.316
beta3_pelagic[7]	29.394	6.617	18.474	29.332	40.333
beta3_pelagic[8]	29.577	6.553	18.602	29.716	40.402
beta3_pelagic[9]	29.481	6.624	18.627	29.562	40.443
beta3_pelagic[10]	29.421	6.700	18.496	29.442	40.391
beta3_pelagic[11]	30.314	1.332	28.147	30.122	33.640
beta3_pelagic[12]	30.000	0.000	30.000	30.000	30.000
beta3_pelagic[13]	30.974	3.763	24.769	30.670	38.563
beta3_pelagic[14]	34.085	1.599	31.470	33.939	37.889
beta3_pelagic[15]	34.127	5.770	23.478	34.954	41.766
beta3_pelagic[16]	30.459	3.023	24.924	30.533	37.200
beta0_yellow[1]	-0.483	0.078	-0.636	-0.484	-0.330
beta0_yellow[2]	0.527	0.099	0.284	0.538	0.677
beta0_yellow[3]	-0.218	0.137	-0.502	-0.218	0.012
beta0_yellow[4]	0.920	0.219	0.387	0.982	1.192
beta0_yellow[5]	0.256	0.328	-0.412	0.271	0.845
beta0_yellow[6]	0.818	0.107	0.617	0.818	1.028
beta0_yellow[7]	0.994	0.091	0.811	0.994	1.170
beta0_yellow[8]	1.025	0.092	0.848	1.022	1.209
beta0_yellow[9]	0.460	0.102	0.262	0.461	0.663
beta0_yellow[10]	0.539	0.071	0.409	0.537	0.679
beta0_yellow[11]	-0.209	0.132	-0.488	-0.197	0.005
beta0_yellow[12]	-1.081	0.095	-1.269	-1.081	-0.902
beta0_yellow[13]	-1.339	0.174	-1.736	-1.323	-1.048
beta0_yellow[14]	-0.252	0.151	-0.588	-0.238	-0.005
beta0_yellow[15]	-1.282	0.259	-1.967	-1.237	-1.027
beta0_yellow[16]	-0.250	0.093	-0.432	-0.250	-0.069
beta1_yellow[1]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[2]	1.038	0.255	0.600	1.022	1.612
beta1_yellow[3]	0.488	0.268	0.021	0.494	0.964
beta1_yellow[4]	1.297	0.545	0.631	1.171	3.066
beta1_yellow[5]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[6]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[7]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[8]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[9]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[10]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[11]	0.332	0.386	0.001	0.172	1.291
beta1_yellow[12]	0.000	0.000	0.000	0.000	0.000
beta1_yellow[13]	1.797	1.631	0.034	1.065	5.837
beta1_yellow[14]	0.703	0.960	0.012	0.431	4.044
beta1_yellow[15]	0.192	0.297	0.011	0.073	1.081
beta1_yellow[16]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[1]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[2]	-5.744	5.590	-20.442	-3.671	-0.264
beta2_yellow[3]	-7.520	5.782	-20.913	-6.174	-0.333
beta2_yellow[4]	-5.843	5.956	-20.771	-3.991	-0.135
beta2_yellow[5]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[6]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[7]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[8]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[9]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[10]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[11]	-7.772	5.828	-21.717	-6.512	-0.433
beta2_yellow[12]	0.000	0.000	0.000	0.000	0.000
beta2_yellow[13]	-5.446	5.200	-18.678	-3.729	-0.289
beta2_yellow[14]	-6.933	5.407	-19.533	-5.706	-0.334
beta2_yellow[15]	-7.896	5.875	-21.184	-6.693	-0.314
beta2_yellow[16]	0.000	0.000	0.000	0.000	0.000
beta3_yellow[1]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[2]	28.985	1.330	27.128	28.750	32.455
beta3_yellow[3]	32.338	3.366	22.635	32.552	38.550
beta3_yellow[4]	29.668	3.177	23.776	28.500	35.538
beta3_yellow[5]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[6]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[7]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[8]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[9]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[10]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[11]	34.749	3.556	30.356	33.768	43.303
beta3_yellow[12]	25.000	0.000	25.000	25.000	25.000
beta3_yellow[13]	32.493	2.665	29.113	32.383	38.039
beta3_yellow[14]	33.922	3.518	29.196	33.528	42.868
beta3_yellow[15]	36.583	3.979	30.376	36.289	43.683
beta3_yellow[16]	35.000	0.000	35.000	35.000	35.000
beta0_black[1]	0.014	0.076	-0.130	0.013	0.167
beta0_black[2]	1.855	0.032	1.792	1.854	1.917
beta0_black[3]	1.239	0.037	1.166	1.238	1.311
beta0_black[4]	2.425	0.048	2.334	2.424	2.520
beta0_black[5]	5.264	2.441	1.910	4.806	11.509
beta0_black[6]	5.314	2.188	2.708	4.804	11.104
beta0_black[7]	4.245	2.266	1.684	3.589	10.425
beta0_black[8]	0.879	0.211	0.488	0.870	1.315
beta0_black[9]	2.818	0.276	2.295	2.812	3.367
beta0_black[10]	1.287	0.051	1.188	1.287	1.388
beta0_black[11]	2.978	0.052	2.877	2.977	3.080
beta0_black[12]	4.914	0.193	4.536	4.913	5.293
beta0_black[13]	0.417	0.318	0.014	0.346	1.200
beta0_black[14]	2.853	0.105	2.649	2.851	3.056
beta0_black[15]	0.723	0.052	0.621	0.723	0.826
beta0_black[16]	3.772	0.074	3.624	3.772	3.910
beta2_black[1]	0.000	0.000	0.000	0.000	0.000
beta2_black[2]	0.000	0.000	0.000	0.000	0.000
beta2_black[3]	0.000	0.000	0.000	0.000	0.000
beta2_black[4]	0.000	0.000	0.000	0.000	0.000
beta2_black[5]	0.000	0.000	0.000	0.000	0.000
beta2_black[6]	0.000	0.000	0.000	0.000	0.000
beta2_black[7]	0.000	0.000	0.000	0.000	0.000
beta2_black[8]	0.000	0.000	0.000	0.000	0.000
beta2_black[9]	0.000	0.000	0.000	0.000	0.000
beta2_black[10]	0.000	0.000	0.000	0.000	0.000
beta2_black[11]	0.000	0.000	0.000	0.000	0.000
beta2_black[12]	0.000	0.000	0.000	0.000	0.000
beta2_black[13]	-2.573	1.841	-7.189	-2.118	-0.330
beta2_black[14]	0.000	0.000	0.000	0.000	0.000
beta2_black[15]	0.000	0.000	0.000	0.000	0.000
beta2_black[16]	0.000	0.000	0.000	0.000	0.000
beta3_black[1]	25.000	0.000	25.000	25.000	25.000
beta3_black[2]	25.000	0.000	25.000	25.000	25.000
beta3_black[3]	25.000	0.000	25.000	25.000	25.000
beta3_black[4]	25.000	0.000	25.000	25.000	25.000
beta3_black[5]	25.000	0.000	25.000	25.000	25.000
beta3_black[6]	25.000	0.000	25.000	25.000	25.000
beta3_black[7]	25.000	0.000	25.000	25.000	25.000
beta3_black[8]	25.000	0.000	25.000	25.000	25.000
beta3_black[9]	25.000	0.000	25.000	25.000	25.000
beta3_black[10]	25.000	0.000	25.000	25.000	25.000
beta3_black[11]	25.000	0.000	25.000	25.000	25.000
beta3_black[12]	25.000	0.000	25.000	25.000	25.000
beta3_black[13]	37.571	2.982	29.448	38.193	40.853
beta3_black[14]	25.000	0.000	25.000	25.000	25.000
beta3_black[15]	25.000	0.000	25.000	25.000	25.000
beta3_black[16]	25.000	0.000	25.000	25.000	25.000
beta4_black[1]	-0.328	0.094	-0.511	-0.326	-0.145
beta4_black[2]	0.256	0.054	0.150	0.256	0.359
beta4_black[3]	-0.719	0.064	-0.844	-0.719	-0.594
beta4_black[4]	0.629	0.128	0.388	0.630	0.885
beta4_black[5]	0.781	1.663	-1.173	0.491	4.373
beta4_black[6]	0.860	1.555	-0.883	0.532	4.653
beta4_black[7]	0.747	1.258	-0.898	0.509	3.912
beta4_black[8]	0.088	0.281	-0.508	0.102	0.606
beta4_black[9]	1.039	0.821	-0.082	0.854	3.067
beta4_black[10]	0.074	0.067	-0.059	0.074	0.204
beta4_black[11]	-0.608	0.065	-0.741	-0.607	-0.485
beta4_black[12]	1.237	0.294	0.689	1.232	1.833
beta4_black[13]	-1.001	0.076	-1.149	-1.000	-0.852
beta4_black[14]	0.224	0.176	-0.113	0.223	0.573
beta4_black[15]	0.002	0.047	-0.090	0.002	0.096
beta4_black[16]	-0.415	0.118	-0.644	-0.414	-0.181
beta0_dsr[11]	-0.329	0.009	-0.347	-0.329	-0.311
beta0_dsr[12]	1.409	0.045	1.324	1.408	1.496
beta0_dsr[13]	0.155	0.021	0.114	0.154	0.197
beta0_dsr[14]	-0.127	0.017	-0.159	-0.127	-0.093
beta0_dsr[15]	0.206	0.014	0.179	0.206	0.233
beta0_dsr[16]	0.159	0.016	0.128	0.159	0.189
beta1_dsr[11]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[12]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[13]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[14]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[15]	0.000	0.000	0.000	0.000	0.000
beta1_dsr[16]	0.000	0.000	0.000	0.000	0.000
beta4_dsr[11]	0.544	0.035	0.472	0.543	0.613
beta4_dsr[12]	0.799	0.227	0.379	0.793	1.246
beta4_dsr[13]	-0.077	0.041	-0.157	-0.078	0.003
beta4_dsr[14]	-0.208	0.070	-0.345	-0.208	-0.074
beta4_dsr[15]	0.217	0.022	0.174	0.217	0.258
beta4_dsr[16]	-0.049	0.046	-0.138	-0.051	0.040
beta0_slope[11]	-2.209	0.021	-2.249	-2.209	-2.167
beta0_slope[12]	-4.985	0.184	-5.368	-4.974	-4.642
beta0_slope[13]	-1.501	0.034	-1.571	-1.500	-1.435
beta0_slope[14]	-2.712	0.051	-2.814	-2.710	-2.613
beta0_slope[15]	-1.592	0.022	-1.637	-1.592	-1.550
beta0_slope[16]	-2.647	0.041	-2.726	-2.648	-2.564
beta1_slope[11]	0.000	0.000	0.000	0.000	0.000
beta1_slope[12]	0.000	0.000	0.000	0.000	0.000
beta1_slope[13]	0.000	0.000	0.000	0.000	0.000
beta1_slope[14]	0.000	0.000	0.000	0.000	0.000
beta1_slope[15]	0.000	0.000	0.000	0.000	0.000
beta1_slope[16]	0.000	0.000	0.000	0.000	0.000
beta4_slope[11]	-0.333	0.076	-0.487	-0.330	-0.186
beta4_slope[12]	-4.667	0.805	-6.337	-4.609	-3.198
beta4_slope[13]	0.326	0.054	0.223	0.326	0.432
beta4_slope[14]	-0.161	0.152	-0.475	-0.156	0.125
beta4_slope[15]	-0.359	0.033	-0.425	-0.359	-0.294
beta4_slope[16]	-0.147	0.102	-0.350	-0.148	0.045
beta5_slope[11]	3.068	0.109	2.866	3.063	3.294
beta5_slope[12]	2.002	0.388	1.287	1.985	2.788
beta5_slope[13]	0.885	0.078	0.731	0.886	1.040
beta5_slope[14]	3.749	0.358	3.109	3.730	4.498
beta5_slope[15]	1.307	0.072	1.167	1.307	1.447
beta5_slope[16]	2.730	0.233	2.290	2.728	3.192
tau_comp[1]	3.770	0.297	3.206	3.760	4.369
tau_comp[2]	2.853	0.250	2.383	2.842	3.387
tau_comp[3]	0.868	0.074	0.727	0.865	1.014
tau_comp[4]	1.786	0.274	1.303	1.765	2.372
tau_comp[5]	1.273	0.186	0.943	1.262	1.676
sigma_H[1]	0.206	0.049	0.121	0.202	0.311
sigma_H[2]	0.176	0.029	0.128	0.174	0.244
sigma_H[3]	0.189	0.041	0.118	0.186	0.278
sigma_H[4]	0.461	0.081	0.323	0.452	0.648
sigma_H[5]	0.959	0.201	0.589	0.948	1.405
sigma_H[6]	0.526	0.194	0.159	0.519	0.931
sigma_H[7]	0.310	0.064	0.211	0.299	0.458
sigma_H[8]	0.418	0.094	0.255	0.409	0.625
sigma_H[9]	0.514	0.119	0.329	0.495	0.791
sigma_H[10]	0.210	0.039	0.142	0.207	0.298
sigma_H[11]	0.267	0.045	0.192	0.262	0.366
sigma_H[12]	0.379	0.147	0.196	0.331	0.712
sigma_H[13]	0.219	0.037	0.158	0.215	0.301
sigma_H[14]	0.495	0.089	0.339	0.489	0.689
sigma_H[15]	0.232	0.038	0.165	0.228	0.317
sigma_H[16]	0.217	0.041	0.147	0.213	0.311
lambda_H[1]	2.301	3.508	0.088	1.056	11.787
lambda_H[2]	5.219	5.195	0.513	3.833	17.868
lambda_H[3]	5.630	8.004	0.254	2.899	28.466
lambda_H[4]	0.007	0.005	0.001	0.005	0.018
lambda_H[5]	2.474	5.454	0.029	0.726	17.861
lambda_H[6]	5.695	12.543	0.008	0.206	43.119
lambda_H[7]	0.012	0.009	0.002	0.010	0.034
lambda_H[8]	7.478	9.456	0.101	4.262	33.928
lambda_H[9]	0.015	0.010	0.003	0.013	0.041
lambda_H[10]	0.226	0.235	0.027	0.156	0.788
lambda_H[11]	0.246	0.342	0.011	0.117	1.188
lambda_H[12]	3.492	4.367	0.161	2.144	15.030
lambda_H[13]	3.225	2.846	0.300	2.426	10.891
lambda_H[14]	3.438	4.209	0.213	2.192	14.623
lambda_H[15]	0.019	0.043	0.003	0.014	0.061
lambda_H[16]	0.610	0.917	0.035	0.301	3.206
mu_lambda_H[1]	4.103	1.875	1.074	3.869	8.167
mu_lambda_H[2]	3.617	1.968	0.478	3.459	7.846
mu_lambda_H[3]	3.386	1.822	0.731	3.091	7.531
sigma_lambda_H[1]	8.137	4.309	1.877	7.346	18.085
sigma_lambda_H[2]	7.987	4.770	0.822	7.319	18.497
sigma_lambda_H[3]	6.150	3.990	0.941	5.235	16.154
beta_H[1,1]	6.543	1.270	3.591	6.767	8.411
beta_H[2,1]	9.758	0.613	8.371	9.796	10.852
beta_H[3,1]	7.822	0.778	6.016	7.924	9.097
beta_H[4,1]	10.472	7.579	-4.911	10.522	24.769
beta_H[5,1]	-0.254	2.531	-5.618	-0.116	4.070
beta_H[6,1]	3.007	4.445	-7.917	4.505	8.473
beta_H[7,1]	1.135	5.683	-10.912	1.548	11.700
beta_H[8,1]	1.257	3.358	-2.577	1.255	3.439
beta_H[9,1]	13.310	5.748	2.282	13.190	25.485
beta_H[10,1]	6.763	1.824	2.843	6.875	10.199
beta_H[11,1]	5.065	3.435	-2.670	5.736	9.909
beta_H[12,1]	2.615	1.100	0.620	2.563	5.012
beta_H[13,1]	9.126	0.862	7.233	9.178	10.612
beta_H[14,1]	2.143	1.045	0.135	2.146	4.268
beta_H[15,1]	-7.040	3.432	-13.241	-7.247	0.171
beta_H[16,1]	3.782	2.890	-0.662	3.375	10.815
beta_H[1,2]	7.784	0.262	7.254	7.791	8.286
beta_H[2,2]	9.882	0.164	9.572	9.881	10.199
beta_H[3,2]	8.796	0.203	8.389	8.798	9.201
beta_H[4,2]	3.068	1.428	0.413	3.043	5.999
beta_H[5,2]	1.778	0.988	-0.215	1.773	3.749
beta_H[6,2]	5.995	1.110	3.550	6.148	7.774
beta_H[7,2]	2.397	1.053	0.457	2.353	4.585
beta_H[8,2]	3.037	1.008	1.385	3.162	4.260
beta_H[9,2]	3.243	1.099	1.082	3.201	5.548
beta_H[10,2]	8.046	0.354	7.330	8.063	8.745
beta_H[11,2]	9.793	0.622	8.840	9.683	11.169
beta_H[12,2]	3.951	0.376	3.261	3.931	4.750
beta_H[13,2]	9.156	0.238	8.713	9.148	9.667
beta_H[14,2]	4.021	0.346	3.364	4.026	4.689
beta_H[15,2]	11.557	0.616	10.275	11.584	12.683
beta_H[16,2]	4.481	0.819	2.998	4.463	6.141
beta_H[1,3]	8.449	0.264	7.987	8.439	8.985
beta_H[2,3]	9.859	0.149	9.565	9.857	10.170
beta_H[3,3]	9.455	0.162	9.126	9.453	9.782
beta_H[4,3]	-2.174	0.865	-3.830	-2.183	-0.429
beta_H[5,3]	3.882	0.595	2.673	3.912	5.022
beta_H[6,3]	8.392	1.181	6.668	8.120	10.878
beta_H[7,3]	-2.558	0.623	-3.803	-2.548	-1.392
beta_H[8,3]	5.307	0.504	4.705	5.229	6.351
beta_H[9,3]	-2.730	0.713	-4.106	-2.743	-1.340
beta_H[10,3]	8.545	0.270	8.013	8.544	9.082
beta_H[11,3]	8.519	0.287	7.904	8.539	9.039
beta_H[12,3]	5.269	0.334	4.491	5.313	5.779
beta_H[13,3]	8.859	0.174	8.501	8.862	9.191
beta_H[14,3]	5.790	0.273	5.187	5.814	6.269
beta_H[15,3]	10.276	0.294	9.708	10.270	10.858
beta_H[16,3]	6.030	0.618	4.742	6.060	7.116
beta_H[1,4]	8.102	0.188	7.711	8.110	8.450
beta_H[2,4]	9.850	0.129	9.584	9.857	10.079
beta_H[3,4]	9.952	0.153	9.626	9.963	10.217
beta_H[4,4]	11.819	0.464	10.865	11.835	12.707
beta_H[5,4]	5.555	0.724	4.349	5.485	7.223
beta_H[6,4]	7.133	0.914	5.247	7.324	8.468
beta_H[7,4]	8.542	0.348	7.821	8.542	9.207
beta_H[8,4]	6.727	0.262	6.219	6.734	7.199
beta_H[9,4]	7.438	0.454	6.518	7.445	8.327
beta_H[10,4]	7.525	0.216	7.117	7.518	7.959
beta_H[11,4]	9.515	0.187	9.147	9.513	9.891
beta_H[12,4]	7.260	0.196	6.886	7.256	7.662
beta_H[13,4]	9.071	0.138	8.797	9.076	9.324
beta_H[14,4]	7.843	0.214	7.443	7.835	8.284
beta_H[15,4]	9.477	0.220	9.050	9.477	9.893
beta_H[16,4]	9.399	0.242	8.948	9.394	9.883
beta_H[1,5]	8.907	0.152	8.590	8.907	9.202
beta_H[2,5]	10.676	0.100	10.483	10.675	10.877
beta_H[3,5]	10.755	0.161	10.460	10.748	11.088
beta_H[4,5]	8.344	0.487	7.406	8.326	9.339
beta_H[5,5]	5.194	0.598	3.672	5.259	6.156
beta_H[6,5]	8.963	0.628	8.014	8.856	10.324
beta_H[7,5]	6.587	0.347	5.926	6.582	7.306
beta_H[8,5]	8.249	0.201	7.919	8.234	8.641
beta_H[9,5]	8.197	0.468	7.251	8.201	9.116
beta_H[10,5]	10.039	0.222	9.604	10.039	10.476
beta_H[11,5]	11.584	0.218	11.159	11.583	12.004
beta_H[12,5]	8.611	0.186	8.223	8.613	8.974
beta_H[13,5]	10.142	0.130	9.886	10.139	10.401
beta_H[14,5]	9.235	0.213	8.831	9.230	9.675
beta_H[15,5]	11.302	0.236	10.833	11.303	11.761
beta_H[16,5]	9.844	0.182	9.463	9.853	10.177
beta_H[1,6]	10.114	0.231	9.736	10.084	10.665
beta_H[2,6]	11.473	0.118	11.244	11.473	11.702
beta_H[3,6]	10.618	0.156	10.287	10.626	10.903
beta_H[4,6]	12.908	0.883	11.142	12.928	14.558
beta_H[5,6]	5.654	0.607	4.592	5.624	6.936
beta_H[6,6]	8.715	0.578	7.423	8.749	9.788
beta_H[7,6]	10.537	0.590	9.341	10.552	11.705
beta_H[8,6]	9.527	0.253	9.048	9.538	9.954
beta_H[9,6]	8.278	0.785	6.738	8.268	9.847
beta_H[10,6]	9.152	0.346	8.418	9.178	9.764
beta_H[11,6]	10.711	0.360	9.947	10.736	11.340
beta_H[12,6]	9.377	0.237	8.939	9.367	9.869
beta_H[13,6]	11.092	0.148	10.805	11.090	11.389
beta_H[14,6]	9.918	0.279	9.338	9.923	10.458
beta_H[15,6]	10.571	0.422	9.701	10.573	11.408
beta_H[16,6]	10.605	0.248	10.100	10.609	11.089
beta_H[1,7]	10.206	1.204	7.164	10.500	11.877
beta_H[2,7]	12.104	0.532	10.919	12.141	13.070
beta_H[3,7]	10.295	0.665	8.854	10.351	11.403
beta_H[4,7]	1.468	4.636	-7.299	1.352	10.630
beta_H[5,7]	5.865	1.908	2.174	5.912	9.702
beta_H[6,7]	7.023	3.213	-1.000	8.116	10.896
beta_H[7,7]	7.694	2.948	2.037	7.668	13.544
beta_H[8,7]	10.863	0.968	9.196	10.859	12.481
beta_H[9,7]	5.070	3.991	-3.007	5.129	12.720
beta_H[10,7]	9.994	1.613	7.228	9.833	13.550
beta_H[11,7]	11.338	1.793	8.160	11.165	15.334
beta_H[12,7]	9.762	0.999	7.596	9.849	11.396
beta_H[13,7]	11.720	0.708	10.052	11.789	12.924
beta_H[14,7]	10.633	0.892	8.799	10.678	12.250
beta_H[15,7]	13.862	2.204	9.554	13.864	18.267
beta_H[16,7]	12.738	1.362	10.353	12.610	15.718
beta0_H[1]	9.149	17.963	-26.037	8.939	45.422
beta0_H[2]	10.524	8.076	-4.773	10.545	27.092
beta0_H[3]	9.631	10.734	-10.702	9.734	30.369
beta0_H[4]	9.475	176.716	-345.977	7.494	364.342
beta0_H[5]	4.721	27.708	-46.460	4.318	57.761
beta0_H[6]	8.192	52.992	-98.212	7.771	129.531
beta0_H[7]	5.820	143.553	-288.824	6.375	308.591
beta0_H[8]	6.196	24.583	-19.467	6.430	31.705
beta0_H[9]	5.741	123.768	-249.177	5.480	251.950
beta0_H[10]	8.333	34.744	-63.328	8.964	76.485
beta0_H[11]	9.531	50.108	-100.927	9.882	114.943
beta0_H[12]	6.522	13.098	-18.760	6.337	32.405
beta0_H[13]	9.781	10.390	-11.381	9.762	29.457
beta0_H[14]	7.102	11.019	-15.804	7.037	27.830
beta0_H[15]	9.147	116.228	-225.346	9.307	252.444
beta0_H[16]	7.239	28.399	-51.772	7.413	64.850